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Background: The reaction ^^0(p, 7)^^F influences hydrogen-burning nucleosynthesis in several stellar sites, such 
as red giants, asymptotic giant branch (AGE) stars, massive stars and classical novae. In the relevant temperature 
range for these environments (Tg = 0.01-0.4), the main contributions to the rate of this reaction are the direct 
capture process, two low lying narrow resonances {Er = 65.1 and 183 keV) and the low-energy tails of two broad 
resonances [Er — 557 and 677 keV). 

Purpose: Previous measurements and calculations give contradictory results for the direct capture contribution 
which in turn increases the uncertainty of the reaction rate. In addition, very few published cross section data 
exist for the high energy region that might affect the interpretation of the direct capture and the contributions 
of the broad resonances in the lower energy range. This work aims to address these issues. 

Method: The reaction cross section was measured in a wide proton energy range yhjc.m. = 345 - 1700 keV) and 
at several angles {Oiah = 0°, 45°, 90°, 135°). The observed primary 7-transitions were used as input in an i?-matrix 
code in order to obtain the contribution of the direct capture and the two broad resonances to the low-energy 
region. 

Results: The extrapolated S-factor from the present data is in good agreement with the existing literature data 
in the low-energy region. A new reaction rate was calculated from the combined results of this work and literature 
S-factor determinations. Resonance strengths and branchings are reported for several ^^F states. 

Conclusions: We were able to extrapolate the astrophysical S-factor of the reaction ^^0(p, 7)^^F at low energies 
from cross section data taken at higher energies. No significant changes in the nucleosynthesis are expected from 
the newly calculated reaction rate. 



I. INTRODUCTION 

Proton induced reactions on ^^O nuclei take place 
in the hydrogen burning shells of red giants, asymp- 
totic giant branch (AGB) stars, and the cores of 
massive stars, at temperatures around T = 0.03 — 0.1 
GK. At these conditions hydrogen burning is domi- 
nated by the CNO cycles [1], namely the reaction se- 
quences i2c(p, 7)i^N(/3+, vY^C{p, 7)^'^N(p, 7)^^0(/3+, v) 
i5N(p,a)i2C (CNO-I), i^N(p,7)i60(p,7)i^F(/3+,^)i^O 
(p, aY^^ (CNO-II), i^O(p, 7)^^F(/3+, vY^O{p, af^^ 
(CNO-III), and ^^0{p,-fY^¥{p,aY^O (CNO-IV). Of 
particular interest is the branching between the two 
reactions "^^0(p, 7)"'^^F and ^^0{p^aY^^ that determines 
the leakage towards the third and fourth CNO cycles. 
The reaction rates of these two proton induced reactions 
determines the branching ratio which in turn affects the 
nucleosynthesis and the abundance ratio of the rare ^^O 
and ^^O isotopes at different environmental conditions 
of the burning site. Comparison with observed abun- 
dance distributions of the oxygen isotopes will provide 
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information on the interplay between nucleosynthesis 
and mixing processes at different stellar burning sites 

The two reactions play a similarly important role in 
the hot- CNO cycles during explosive hydrogen burning 
in classical novae. At these conditions hydrogen burn- 
ing reaches peak temperatures around T= 0.1 — 0.4 GK 
[7 . In this scenario, ^^0(p, 7)^^F directly affects the pro- 
duction of the /3+ unstable ^^F (ti/2 = 110 min) [8, 9j, 
whose 511 keV electron-positron annihilation 7-ray could 
potentially be detected by 7-ray satellites, such as the 
INTEGRAL observatory [10]. 

The stellar reaction rates are determined by the reac- 
tion cross section at the stellar energy range. At low ener- 
gies, the reaction rate of ^^0(p,7)^^F {Q = 5.606 MeV) 
is affected by two low lying narrow resonances (F < 1 
keV) at Er = 65.1 keV and Er = 183 keV {Ep = 68.9 
and 193 keV respectively), which have been subject to 
a number of recent low-energy studies. The resonance 
at Ep = 68.9 keV is extremely hard to measure directly 
with current techniques, and its strength is estimated in- 
directly by experimental and theoretical constrains on its 
partial widths. On the other hand, recent work by Fox et 
al [11 and Chafa et al. [12] have successfully measured 
directly the strength of the resonance at Ep = 193 keV. 

This paper reports on a new measurement of the 
^^0(p,7)^^F reaction. The measurement focused on ex- 
amining the nature and strength of cross section compo- 
nents which extend into the low-energy range and con- 



2 



tribute substantially to the stellar rate. These contribu- 
tions include a strong direct capture component and the 
low-energy tails of two broad resonances at = 557 
keV and Er = 677 keV {Ep = 590 and 717 keV, respec- 
tively). These components had first been investigated by 
Rolfs [13]. The more recent studies by Fox et al [TT] 
and Chafa et al. \12i contradict the earlier cross section 
data and differ by up to a factor of 2. Measurements by 
Newton et al [M] and Hager et al [15] focused on the 
direct capture component at low energies, in the range 
Ec.m. = 250 — 500 keV. The limited energy range in these 
studies prohibits a direct normalization to the strength 
of the higher energy resonances. Newton et al. esti- 
mated the contribution of the two broader resonances 
in their measured energy range by multiplying the pre- 
viously recorded resonance strengths [13 by a factor of 
0.62, following an earlier suggestion by Fox et al. pT] . 

The inconsistencies between the different data sets im- 
pacts the extrapolation of the S'-factor to the stellar en- 
ergy range, which makes a further re-investigation of 
the reaction necessary. The present study aims at con- 
straining the low-energy contributions of the high en- 
ergy resonances and the direct capture by measuring 
the capture cross section of individual primary transi- 
tions over a wide range of energies, Ep = 365 — 1800 
keV {Ec.m. = 345 - 1700 keV) and angles, Oiab = 0°, 
45°, 90°, and 135°. The data are then fitted within the 
framework of a multi-level, multi-channel i?-matrix ap- 
proach [16 , allowing for a more consistent and better- 
constrained extrapolation to lower energies. The angular 
distribution information is important in constraining the 
direct capture contributions of different final state orbital 
angular momenta, when more than one are possible, as 
well as the magnitude of the background poles, when 
these were included. Target yield deconvolution and an- 
gular attenuation effects are carefully taken into account 
and coincidence summing corrections are kept to a min- 
imum, in order to obtain more reliable data. Strengths 
and branching ratios are provided for several resonances 
in this energy region, and are compared to literature val- 
ues. 



Sec. [TTj and III describe the experimental equipment, 
procedure and results. In Sec. [TV] we discuss the R- 



matrix fit of the data and the extrapolation towards lower 
energies in comparison with previous studies. The new 
reaction rate is presented in Sec. [Vj Throughout the pa- 
per, Ecm. refers to the reaction energy at the center-of- 
mass, Ep refers to the proton beam energy in the labora- 
tory reference frame, E^ is the excitation energy of the 
nucleus and Er is the resonance energy in the center-of- 
mass, unless otherwise noted. 



II. EXPERIMENTAL SETUP AND 
PROCEDURE 

The experiments were carried out at the Nuclear Sci- 
ence Laboratory, at the University of Notre Dame. Pro- 
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FIG. 1. Mechanical drawing of the target chamber used for 
the experiments. 



ton beams were provided by the 4 MV KN Van de Graaff 
accelerator covering the energy range Ep = 600 — 1800 
keV. For the study of the lower energy range, Ep = 
365-700 keV the 1 MV JN Van de Graaff accelerator was 
used. Typical beam currents between 20 and 40 /iA were 
achieved for both machines, with a beam energy resolu- 
tion of approximately 1.5 keV and an energy uncertainty 
of 0.5 keV. The proton beam entered the target chamber 
through a liquid-nitrogen-cooled copper pipe (i.e. cold 
trap) in order to reduce carbon build-up on the target 
surface (Fig. [T]). The cold trap was biased to —400 V to 
suppress secondary electrons from the target. The target 
and chamber formed a Faraday cup for charge integra- 
tion. The position of the beam on the target was defined 
by a set of vertical and horizontal slits. The beam was 
swept horizontally and vertically across a target area of 
2 cm^ by steerers to dissipate power over a large target 
area. The target was directly water cooled using deion- 
ized water. 

Two ^^O targets were prepared by anodization of 0.5 
mm sheets of tantalum in 90.1% ^^0-enriched wateij^ 
targets Ta205 "Thin" and "Thick". The former was used 
in the proton energy range Ep = 500 — 1150 keV, and 
the latter between Ep = 365 — 1625 keV. Such targets 
have been shown [22] to have a well-known stoichiom- 
etry (Ta2 05) and their thickness can be controlled by 
the anodizing voltage. In addition, an implanted target 
was prepared by bombarding a Ta backing with a 30 keV 
^^O beam. The implanted target was used in the proton 
energy range Ep = 400 — 180 keV. Several scans of nar- 
row resonances were performed in order to measure the 
thickness of all the targets (Fig.[2|. The thickness of tar- 
get Ta2 05 "Thick" , prepared with the higher anodization 
voltage, was determined by the width of the 7-ray yield 
curves taken for the narrow resonances at Ep = 519 keV 
{AEp = 12.5 ± 0.8 keV) and Ep = 1098 keV (A^;^ = 



1 According to the supplier 9.5% ^^O, 90.1% ^'^O, and 0.4% ^^O. 
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FIG. 2. (Color online) Yield as a function of energy of the 
2100 keV 7-ray from the Ep = 1098 keV {Er = 1037 keV) 
"^^0(p, 7)"^^F narrow resonance. The three sets of data points 
correspond to the three different targets used in the experi- 
ments. 



8.45 ± 0.55 keV). Using stopping power tables [23] and 
the known stoichiometry, the target thickness of Ta205 
"Thick" was measured to be 6.38±0.37x 10^^ atoms/cm^. 
The uncertainty of the thick target was calculated from 
the uncertainty of the measured energy loss (~ 5%) 
and the tabulated stopping powers (~ 4%). The un- 
certainty in the ^^O water enrichment was also taken 
into account (~ 3%). The thicknesses of the other an- 
odized target, Ta205 "Thin", and the implanted tar- 
get were calculated relative to the thickness of Ta205 
"Thick", with values of 3.00 ± 0.32 x 10^^ atoms/cm^ 
and 2.25 ±0.26 x 10^^ atoms/cm^, respectively. The un- 
certainties for these targets arise mainly from the un- 
certainty of the thick target 6%) and the uncertain- 
ties of the measured integrated yields of the resonance 
scans (~ 8%). Resonance scans were repeated periodi- 
cally throughout the experiments to monitor the quality 
of the targets. No target degradation was observed. Due 
to the unknown stoichiometry of the implanted target, 
only data from the high energy region {Ep > 1250 keV) 
where the cross section varies smoothly with beam en- 
ergy were considered in the analysis. The use of the two 
anodized targets with different thicknesses helped to con- 
firm the yield deconvolution method, as discussed later. 

At this point, it is important to note some of the lim- 
iting factors in the detection of 7-rays in the present 
experiment. First, the decay scheme of the compound 
nucleus, ^^F, is rather fragmented, which means that an 
appreciable percentage of decay paths consist of multiple 
transitions with low branching ratios. As a result, some 
transitions may have low statistics and may be unde- 
tectable because of the large Compton continuum from 
the stronger transitions. The latter is even more pro- 
nounced when using crystals of small volume, as in this 
work. Therefore, obtaining a full excitation curve of the 



two broader resonances, Ep = 590 and 717 keV, including 
their tails, for all the branchings was not possible in the 
timeframe of the experiments, and given the background 
involved. Instead, long runs on-resonance were used to 
obtain the branching ratio information, as discussed in 
the next section. 

At 7-ray energies Ej = 5.0 — 6.2 MeV, the spectrum 
is dominated by the presence of the 6130 keV line from 
the ^^F(p, a7)^^0 reaction, whose cross section exhibits 
several strong resonances. This is a common problem 
when using tantalum backings, since fluorine is part of 
the extraction process of tantalum, making it very hard 
to obtain tantalum with sufficiently low fluorine contam- 
ination levels. Heating of the tantalum backing in a high 
vacuum chamber did not decrease the ^^F contamination 
levels noticeably. For this reason, the detection of the 
primaries R/DC^ 937 and R/DC^ 1121 was inhibited 
for a wide proton energy range but more importantly at 
energies Ep = 1250 - 1450 keV {Ec.m. = 1180 - 1370 
keV), where the ^^F{p^ aj^^O reaction exhibits strong 
resonances. ^^O contamination did not affect the mea- 
surements significantly, since the enriched water used for 
the production of the anodized targets contained only 
a very small percentage of ^^O. Cross section plots are 



presented in Sec. [IV|and in the appendix along with the 
i?-matrix fits. 

Angular distributions of the prompt 7-rays from the 
^^0(p, 7)^^F reaction were measured with four 20% (rel- 
ative to a 3" X 3" Nal detector) high-purity Ge (HPGe) 
detectors placed at 0°, 45°, 90°, 135° at 8 cm distance 
from the center of the target (Fig.[3|i). Two typical spec- 
tra showing the primary transition to the first excited 
state in ^^F {E^, = 937 keV) obained with the 45° de- 
tector are shown in Fig. [4j The spectrum in Fig. [4^ was 
recorded at proton energy Ep = 1625 keV {Ec.m. = 1535 
keV), where the cross section is dominated by the direct 
capture mechanism with an approximate value of ~ 200 
nb/sr. The spectrum in Fig. |4]3 was recorded at proton 
energy Ep = 400 keV {Ec.m. — 378 keV), which corre- 
sponds to the low-energy side of the Ep = 590 keV res- 
onance with an approximate cross section of ~ 4 nb/sr. 
In addition, measurements were performed with a HPGe 
clover detector, segmented in four 25% crystals, placed 




FIG. 3. Detector geometries of 
setup and b) the clover setup. 



the angular distribution 
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FIG. 4. Spectra of the 45° detector of the angular distribution setup in the energy region of the primary transition to the first 
excited state in ^^F (Ea^ = 937 keV) for proton beam energies of a) 1625 keV (Ec.m. = 1535 keV) and b) 400 keV {Ec.m. = 378 
keV). The low counting statistics of the low-energy spectrum is indicative of the low cross sections, in combination with 
the detector geometry. The more prominent peaks in both spectra are due to the ^^F(p, a^y^O reaction, from the fluorine 
contamination in the tantalum backing. 



at 45° and 1 cm from the target (Fig.jsjD). The spectra of 
the four individual crystals were added off-line and were 
analyzed as one. Performing the experiment with two dif- 
ferent detector geometries was important in testing the 
validity of the data analysis procedure. The intrinsic en- 
ergy resolution for all crystals is 2.2 keV at = 1.33 
MeV. 

Full peak detection efficiencies as functions of the 7- 
ray energy were obtained for both setups using cal- 
ibrated (^^Co, i^^Cs) and uncalibrated (^^Co) 7-ray 
sources as well as the well-known narrow resonances 
of 27Al(p,7)28si and ^^N{p,jy^O at Ep = 992 keV 
{Ec.m. = 957 keV) [17J and 278 keV {Ec.m. = 260 
keV) [18], respectively. Detection efficiencies for the sin- 
gle and double escape peaks were also measured, since 
most primaries of interest exhibited strong escape lines 
due to the high Q- value of the reaction. Figure |5] shows 
an example of a full energy peak, a single escape peak 
and a double escape peak efficiency curve. The peak ef- 
ficiency curve was fitted with the function [19] 

r]pe = exp(a + b ln{E^) + c ln(£:^)^ + d ln{E^)^), (1) 

whereas for the single and double escape curves sim- 
ple second order polynomials were used. Coincidence 
summing-in corrections were negligible, given the low 
peak efficiencies (see Fig. |5|. A clear indication of this 
was the fact that no peak appeared in the aluminum res- 
onance runs in the 7-ray spectra of the clover detector 
setup at the 7-energy of 12541 keV, where a summing-in 
peak usually appears from the sum of the strong 1779 
and 10761 keV transitions in ^^Si. 

Coincidence summing-out corrections for each primary 
transition were calculated using a Monte Carlo program 
that required input for the total efficiency of the detec- 
tor set-up as a function of energy and the level scheme of 
-'^^F (taken from [261 [27]). The program would simulate 



a million decays and subtract from the 7-ray of interest 
according to the coincidence 7-rays and the total effi- 
ciency information. The total efficiency was measured at 
two energy points with the ^^Co and ^^^Cs sources and 
was extrapolated to the entire energy range of interest 
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FIG. 5. (Color online) Peak efficiency curves for the 0° de- 
tector from the angular distribution experiment. The red 
circles refer to 7 strength measurements of transitions in 
the ^^Al(p, 7)^^Si reaction, the blue squares refer to 7 decay 
measurements with ^°Co and ^^^Cs sources, the black trian- 
gles relate to 7 strength measurements associated with the 
^'^N(p, 7)^^0 reaction and the open red circles refer to ^^Co 
source measurements. The two dashed lines correspond to the 
single (purple dot-dashed) and double (green dashed) escape 
peak efficiencies. The bottom graph shows the deviation of 
the full peak efficiency fit from the data points in units of 
standard deviation. 
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using the Monte Carlo simulation package Geant4 [20], 
normalized to the two data points. Typical total effi- 
ciency values for the angular distribution setup were ap- 
proximately 1.5% and for the clover setup 5%. Due to 
the low total efficiencies, summing-out corrections were 
small, and their contribution to the cross section uncer- 
tainty was < 1%. 

Angular attenuation coefficients were calculated from 
the formula |21] 
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Qn 



^ J^-'' Pnjcos P)7i{l3,E) sin 13 dl3 
Jl^--T1{I3,E) sin 13 dl3 



(2) 



with /3 being the angle relative to the detector's symme- 
try axis and iri{j3^E) the peak efficiency as a function of 
angle and energy ^ . Qn does not depend strongly on 
the absolute scale of the peak efficiency and as a conse- 
quence on the energy of the 7-ray. For the present de- 
tector setups, 7^(/3, E) was estimated from Geant4 sim- 
ulations. The results were Qi = 0.988, Q2 = 0.963, 
Qs = 0.928 and (^4 = 0.880 for the angular distribution 
setup and Qi = 0.755, Q2 = 0.388, Q3 = 0.067 and 
Q4 = —0.103 for the clover setup. 



III. ANALYSIS AND RESULTS 

All detectable primary transitions from the 
""^^©(p, 7)^^F reaction and the secondary transition 
from the ffist excited state of ^^F at 937 keV to the 
ground state were measured in the proton energy range 
from 365 to 1800 keV. The data from the primaries were 
used in the i?-matrix analysis. Secondary transitions 
were not included in the i?-matrix analysis but were im- 
portant for the direct comparison to the Rolfs data. The 
experimental differential yield for any 7-ray transition 
was calculated with the equation 



(3) 



where Ep is the energy of the proton beam in the labo- 
ratory system, Oiab is the angle of the detector, Ne^ is 
the number of counts of the 7-ray peak and Nproj are the 
number of protons incident on target. The last quantity 
is calculated from the total charge collected on the target 
throughout each run. The deconvolution of target effects 
from the experimental yield curves was performed by di- 
rectly performing an i?-matrix fit on the experimental 
yields rather than the cross sections, taking into account 
the thickness of the targets and the effective stopping 
power. As a ffist step, each experimental yield curve is 
fitted separately, using reasonable i^-matrix parameter 
input. The experimental cross section is then estimated 
for each measured energy point from the fitted yield curve 
and the corresponding calculated cross section by 
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FIG. 6. (Color online) a) Measured yield of the secondary 
transition from the 1^* excited state to the ground state of ^^F 
(937^G.S.) in the energy region of the two broad resonances. 
The two data sets correspond to the two anodized targets, 
Ta205 "Thin" (red circles) and "Thick" (black triangles), b) 
Converted yield to cross section with the method described 
in the text. Narrow resonances were removed in this analysis. 



In this equation, Yexp{Ei) is the experimental yield given 
in Eq. [sj Ycaic{Ei) is the yield obtained by the i?-matrix 
fit taking into account the target thickness, and (Jcaic{Ei) 
is the cross section that corresponds to the calculated 
yield when the target effects are removed. This method 
cannot be applied for narrow resonances, as these have 
features that cannot be perfectly reproduced by the code, 
such as beam energy resolution and straggling. For this 
reason, the properties of the narrow resonances were 
studied separately. 

Evidence for the validity of the deconvolution method 
is presented in Fig.|6] Fig. [6^ shows the measured yield of 
the secondary transition from the 1st excited state to the 
ground state of ^^F (937^G.S.). The black triangles and 
red circles correspond to the "Thick" and "Thin" target 
yield excitation functions, respectively. Fig. (6)3 shows the 
cross section calculated from the yield with the method 
described above. The convergence of the two curves is 
evidence that the procedure works well. A small sys- 
tematic deviation can be explained by the uncertainties 
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of the target thicknesses. The yields of ah observed pri- 
mary transitions were converted to cross section with this 
method, excluding most narrow resonances. 

The strengths of the narrow resonances at Ep = 519, 
826, 930, and 1098 keV {Er = 490, 780, 878, and 
1037 keV, respectively) were measured with the Ta205 
"Thick" target, to take advantage of its well-known sto- 
ichiometry and thick target yield. Long runs off res- 
onance were also taken in order to identify contribu- 
tions that do not correspond to the measured resonance 
and subtracted from the yield data. The angle inte- 
grated yield of each transition was obtained by fitting 
the differential yields of each detector with the equation 
W{0) = Ymaxi'^ + a2P2{cos6))/47r. The angle integrated 
yields for each of the observed transitions were added to 
obtain the branching ratios and the resonance strengths. 
The latter quantity was calculated using the formula 



^eff,r 



arctan(Ai;/r)' 



(5) 



where Y^ax is the total angle integrated 7-ray yield mea- 
sured on the target yield plateau, is the de Broglie 
wavelength at the resonance energy, is the energy 
loss of the proton beam inside the target at the resonance 
energy, T is the total width of the resonance, taken from 
Ref. [25, and eeff,r is the effective center-of-mass stop- 
ping power of protons in tantalum pentoxide at the reso- 
nance energy [23 . The strengths of the two broader res- 
onances at Ep = 590 keV and 717 keV were determined 



from the i?- matrix fit parameters, presented in Sec. IV 



Tables [l| and [ll] summarize the results for the strengths 
and branching ratios of the six measured resonances. In 
Table [l| the results are compared with literature data. 
The values of Rolfs, corrected by a factor of 0.62 agree 
very well with the present data. This result confirms the 
analysis by Fox et al. [11 indicating that the strengths 
quoted by Rolfs [13] have to be renormalized using the 
updated value for the Ep = 632 keV resonance strength 
in the ^^Al(j9, 7)^^81 reaction [25 . The measurements by 
Sens et al. [27] are systematically lower by a factor of 2 
compared to the present work, possibly due to the use 
of unpublished stopping power values. Finally, it is im- 
portant to note that the resonance strength of the 519 
keV resonance is in excellent agreement with the more 
recent measurement by Fox et al. The experimental un- 
certainties arise primarily from statistics, the uncertainty 
in the absolute (8%) and relative (3%) peak efficiency and 
the uncertainty of the effective stopping power (~ 4%). 
The branching ratios for the strong transitions reported 
in Table [ll] agree very well with values found in litera- 
ture [26] |27] . For weaker transitions some discrepancies 
are observed which could point to unaccounted summing 
effects in the literature values. 

The cross section of the secondary transition from the 
first excited state of -^^F to its ground state (937 ^G.S.) 
was measured to be ~ 40% lower than that of Rolfs (see 
figure 11 of [13 ). Note that this 40% discrepancy is 
not related to the 0.62 factor with which the resonance 



strengths were renormalized, at least to the best of the 
authors' knowledge. It is not clear if or how much of 
this discrepancy is related to the erroneous resonance 
strength value used for the efficiency calibration, since 
the yield of this transition was measured relative to that 
of ^^0(p, 7). To check the self-consistency of the present 
analysis, the yield of this transition was also calculated 
using the measured yield for the primary transitions and 
the known branching through the 937 keV ^^F state. The 
calculated expected yield and the measured yield of the 
937 keV transition were then compared and found to 
agree to within 5% of each other, a strong indication 
of a correct efficiency calibration and branching ratio in- 
formation. 



IV. R-MATRIX ANALYSIS AND DISCUSSION 

A multi-channel, multi-level i?-matrix analysis, as de- 
scribed earlier [16], was used to simultaneously fit all the 
measured primary transitions of the ^^0(p, 7) reaction, 
as well as the ^^0(p, a) reaction measured by Kieser et 
28 . As explained in Sec. 



al. 
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the yield data were 



converted to cross sections by a target effect deconvolu- 
tion method. The results from that analysis were used 
as input for the global i?-matrix fit, as is discussed in 
this section. Data from most narrow resonances were ex- 
cluded from the fits since the required number of target 
integration steps would have slowed down the calculation 
time substantially, without improving the accuracy of the 
fit and extrapolation. The direct capture was assumed to 
be El only, as was suggested by previous studies [11 ,[13]. 
The i?-matrix radius for the proton channel was taken as 
Vc = rox {A^^ + ^y^) = 4.46 fm, with ro = 1.25 fm. 
The i?-matrix radius for the a-channel was taken to be 
ac = 5.44 fm. The choice of ro = 1.25 fm is common 
in literature (e.g. [25]), and was adopted as a reasonable 
estimate. The effect of the channel radius value on the 
results is discussed later. 

In i?-matrix theory the direct capture (DC) part of the 
cross section is divided into an external capture part and 
an internal part. The external capture component has a 
given energy dependence, whereas its magnitude is de- 
fined by the asymptotic normalization coefficient (ANC) 
of the bound state. The ANCs can either be treated as 
parameters of the fit, or they can be derived from exper- 
imental spectroscopic factors {C'^S) of the bound states 
taken from literature. The internal part is treated with 
high energy background poles. 

For each of the observed primary transitions the ANC 
of the final state is one of the i?-matrix parameters. Dif- 
ferent than resonance parameters these parameters are 
only restricted by the respective cross sections. Because 
ANC values are only poorly restricted from the present 
data set, literature values derived from transfer reactions 
[29[ [30] were taken where appropriate. 

To test if these values are consistent with our data, two 
i?-matrix fits were performed: a test fit and a main fit. 
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TABLE I. Resonance strengths in ^'^0(p,7)^^F 


Ep (Er) 




Present 




Literature 




(keV) 


J" 


uj (eV) 


Rolfs [I3] 


Rolfs, scale(|j 


Sens et al 27 


Fox et al [TT\ 


519 (490) 


4~ 


0.0130 ± 0.0015 


0.021 ± 0.004 


0.0130 ±0.0025 


0.0068 ± 0.0020 


0.0137 ± 0.0022 


590 (557) 


3+ 


0.37 ±0.05 


0.56 ±0.15 


0.35 ±0.09 


0.18 ±0.05 




717 (677) 


2+ 


0.58 ±0.07 


0.76 ±0.19 


0.47 ±0.12 


0.206 ±0.065 




826 (780) 


2+ 


0.0323 ± 0.0025 


0.050 ±0.015 


0.0309 ± 0.0093 


0.0175 ± 0.0054 




930 (878) 


3+ 


0.0194 ±0.0018 


0.030 ±0.012 


0.0186 ±0.0074 


0.010 ±0.005 




1098 (1037) 


2- 


0.297 ±0.033 


0.36 ±0.10 


0.22 ±0.07 


0.130 ±0.043 





^ The resonance strengths given by Rolfs have been corrected by a factor of 0.62 as suggested by Fox et al. in |11| . 



TABLE IL Branchings of certain resonances in ^^F 





Ex = 6096 


6163 


6283 


6385 


6485 


6644 


Transition 


Ep = 519 


590 


717 


826 


930 


1098 


R -^0 






38 ± 02 


1.4 ± 0.4 


15.1 ± 0.7 




jl ^937 


4.2 ± 0.5 


50.3 ± 0.6 


65.73 ± 0.35 


74.3 ± 4.3 


28.7 ± 2.0 


6.4 ± 0.5 


R ^1040 






0.69 ±0.07 








i^^ll21 


58.6 ±2.3 


2.1 ±0.2 






14.2 ± 1.7 




R ^1700 






5.67 ±0.13 


8.7 ±0.5 


7.9 ±0.6 




R ^2100 


25.1 ± 1.1 




1.74 ±0.09 






58.4 ±2.7 


R ^2523 




7.0 ±0.2 


0.73 ± 0.06 




4.7 ±0.7 




R ^3062 




2.3 ±0.1 






15.8 ± 1.6 




R ^3133 






1.29 ±0.06 






22.1 ± 1.4 


R ^3358 






2.10 ±0.05 








R ^3724 






2.51 ±0.05 






1.6 ±0.3 


R ^3791 


1.3 ±0.2 


10.2 ±0.2 






4.1 ±0.3 


2.8 ±0.4 


R ^3839 




22.0 ±0.3 


13.28 ±0.10 


11.8 ±0.8 


8.0 ±1.1 




i?^4115 


1.8 ±0.3 


2.2 ±0.3 


3.73 ± 0.06 


3.8 ±0.5 




2.2 ±0.3 


R ^4226 




1.4 ±0.2 










R ^4360 






2.14 ±0.06 








R ^4398 


2.2 ±0.3 


2.5 ±0.1 










R ^4652 


6.8 ±0.3 












R ^4860 












2.2 ±0.3 


R ^4964 










1.6 ±0.2 




R ^5502 












4.3 ±0.3 



In the test fit, no background poles were included and 
the ANCs were left as free parameters. The resulting 
ANCs were converted to spectroscopic factors using the 
equation 

C'S={^. ,6, 

where C'^S is the proton single particle spectroscopic fac- 
tor of a given bound state of ^^F with the isospin Clebsch- 
Gordan coefficient factored in, Cfu is the ANC value ob- 
tained by the fit, and hsp is the ANC value obtained by 
a simple single particle mean field calculation, using a 

1 /3 

Wood-Saxon potential with radius r = 1.25^4^ fm, dif- 
fuseness a = 0.65 fm and depth chosen to reproduce the 



binding energy of the state. The spectroscopic factors 
C'^S were then compared to literature spectroscopic fac- 
tors from trans fer r eactions [29, SOj^to determine their 
validity (Table [ml). If a fitted ANC of a given transi- 
tion corresponds to a spectroscopic factor that is much 
larger than the one found in literature, that indicates 
that the transition cannot be fully described by the ex- 
ternal capture component of the i?-matrix theory, and a 



^ The bound state potential used in the study by ^29^ is identical 
to the one used in this study to calculate the single particle 
spectroscopic factor. The potential used by |30j is similar to 
this study's with the exception of a spin-orbit term. 
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R/DC^ 937 keV 

I . . r I 



0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 

E,„. (MeV) 




(MeV) 



FIG. 7. (Color online) i?-matrix fit of the differential cross 
section of the ^^0(p, 7)^^F*(937 keV). The different colors 
and shapes correspond to different angles. 



background pole is required. 

In the main fit, the spectroscopic factors that disagreed 
with the literature data were fixed close to experimental 
spectroscopic factor values, and the missing contributions 
were replaced by high energy background poles. These 
poles were placed at an excitation energy of 15 MeV, 
which is about 8 MeV above the energy region of the 
experimental data, with spins and parities chosen so that 
primary transitions of interest decay via El. As there are 
no available proton scattering data or other appropriate 
experimental data to constrain the proton partial widths 
of the poles, they were fixed at Tp = 6 MeV, chosen 
close to the Wigner limit given by the formula Tw = 
2Pi{Er)jw^ where Pi{Er) is the penetrability factor for 
a given partial wave at the resonance energy and = 
{3/2)h^/{/irl). The T^ partial width was left as a fit 
parameter. 

The fits of the two strongest transitions, R/DC^ 937 
and R/DC^ 3839 keV are shown in Figs. [7| and [sj Their 
contribution to the low-energy total S'-factor amounts to 
(~ 40%) and (~ 16%) respectively. The fits of the weaker 
transitions are shown in the appendix. The reduced-x^ 
of all the fits ranged between 1.5 and 2.7, with the excep- 
tion of the R/DC^ 937 fit that had a value of 6.5. The 
high reduced-x^ values are due to the fact that this is a 
simultaneous fit of multiple transitions, multiple angles 
and two targets of different thicknesses. In addition, im- 
perfect target effect corrections have a big effect on the 
reduced-x^ when regions of sharp changes in the cross 
section are included, as are in the R/DC^ 937 transi- 
tion. 

A good fit was also obtained in this global fit for the 
''^^0(p, a)^^N, shown in Fig. 14 in the appendix. No tar- 



FIG. 8. (Color online) i?-matrix fit of the differential cross 
section of the ^'^0(p, 7)^^F*(3839 keV). The different colors 
and shapes correspond to different angles. The 45° and 135° 
lines lie on top of each other and cannot be distinguished. 



does not affect the results for the two broad resonances 
at Er = 557 and 677 keV. This was verified by a separate 
fit of the (p, a) data alone with target effects included. 
This led to a better description of the narrow resonances, 
but it did not have an influence on the parameters for the 
broad resonances. 

The fitted and fixed ANC values from all the fits and 
the respective calculated spectroscopic factors are shown 
in Table [nil In bold font are the ANCs that were fixed ac- 
cording to literature spectroscopic factors from [29] . The 
latter values were preferred over those from [30] since the 
nuclear potential used in that work was identical to the 
one used in this work for the conversion of the ANCs to 
spectroscopic factors. The ANC values in normal font 
were treated as free parameters, since the corresponding 
spectroscopic factors were close to the literature values. 
It is clear that the validity of this approach to estimate 
the direct capture component depends on the reliabil- 
ity of the ANC values. New measurements to verify the 
ANC values used in this work would be desirable. In 
the present study, an attempt is made to quantify the 
sensitivity of the extrapolation to the ANC values, as 
discussed below. 

The resulting resonance parameters of the main fit are 
shown in Table 



IV 



get effects were taken into account in the case of the (p, a) 
channel, as that would slow down the calculations. This 
leads to a poor fit for the narrow resonances, however, it 



Values in parentheses denote fixed 
parameters. The energies of the two broad resonances 
are in excellent agreement with literature values. The 
proton width of the 6.16 MeV state, F^ = 14.1 ± 0.3 
keV, agrees very well with the reported literature values 
of 13 ± 1 keV ^ and 14.7 ± 1.5 keV [27 , whereas some 
disagreement is observed for the proton width of the 6.28 
MeV state, Tp = 11.3±0.2 keV, compared to the litera- 
ture values, 8.0 ±2.0 keV [26] and 8.5 ±1.0 keV [27]. The 
reported uncertainties include only the statistical uncer- 
tainties that arise from the fit. In addition, there are sys- 
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TABLE III. ANCs and spectroscopic factors 




Level 




Fixed/Fitte 


J— 1 


Literature C^S 




Energy 


h 


ANC 




m 


|30] 




937 





6.1 


0.32 


0.30 <x < 0.34 


O.IC 


b 


937 


2 


1.2 


0.18 


< 0.36 


0.3C 


b 


1121 


2 


2.7 


1.0 


1.00 


0.8( 




2523 





1.4 


0.025 


0.025 


0.014 


2523 


2 


- 


- 


< 0.005 


0.011 


3062 





4.5 


0.32 


0.16 


0.21 


3062 


2 


1.0 


0.37 


0.74 


0.62 


3839 





4.6 


0.42 


0.50 






3839 


2 


0.6 


0.19 


< 0.56 






4115 





2.5 


0.13 


0.13 


0.17 


4115 


2 


1.0 


0.75 


0.75 


0.68 


4652 


2 


1.3 


1.6 


0.98 


1.04 


4964 





3.2 


0.17 


0.22 


0.17 


4964 


2 


0.7 


0.5 


0.49 


0.51 



^ Values in bold font were adopted from literature (for details see 

discussion in text). 
^ According to [30] these spectroscopic factors have a large 

uncertainty due to a broad minimum. 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 'I-r- l 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 

Ee.^.(MeV) 

FIG. 9. 5'-factor obtained from the i?-matrix fit. The con- 
tinuous line represents the total 5'-factor from the best R- 
matrix fit, with background poles included and some ANC 
values fixed. The dashed-dotted line and the dashed line are 
the contributions of the direct capture (background poles in- 
cluded) and the high energy resonances respectively. Some 
narrow resonances are omitted from this plot. 



tematic errors arising from the choice of fixed i?-matrix 
parameters (e.g. channel radii, ANC). Since these param- 
eters are correlated, these systematic errors are difficult 
to extract. However, deviations greater than 0.2 keV 
were never observed from the various fits. Additional 
background poles were also tested, in order to improve 
the fits of the two broader resonances, but no significant 
effects were observed. 

In order to make the comparison with previous work 
easier, the extrapolated cross section curve was converted 
to the astrophysical S'-factor curve given by the formula 

S{E) = a{E)Ee^^'^, (7) 

1/2 

where r] = 0.1575 x Z1Z2 {/J^/E) ' is the Sommerfeld 
parameter, with \i in units of u, and E in MeV. 

Figure [9] shows the total S'-factor of the reaction ob- 
tained from the main i?- matrix fit (continuous line). The 
figure also shows the resonant and direct capture con- 
tributions separately in dashed and dashed-dotted lines 
respectively. As pointed out by [TT] and [12] , the total di- 
rect capture contribution dominates at energies less than 
Ec.m. = 400 keV, when compared with the total con- 
tribution of the tails of all the high energy resonances. 
No significant interference was observed between the two 
components at the low energy side of the two broad reso- 
nances, as the two mechanisms involve different initial or- 
bital angular momenta {li = 1,3 for direct, k = for the 
resonances, known from angular momentum and parity 
conservation considerations). The two broad resonances 
decay mostly via Ml transitions as opposed to El tran- 
sitions for the direct capture. 

The result of the i?-matrix extrapolation had little de- 



pendence on the choice of the channel radius (~ 3% for 
the R/DC^937 extrapolation for a radius 4 < Tc < 5 
fm) or the position of the background poles (~ 5% for 
the total extrapolation when the pole was placed at 35 
MeV instead of 15 MeV). 




0.15 0.20 0.25 0.30 035 0.40 0.45 

Ee.^. (MeV) 

FIG. 10. (Color online) Comparison of the 5'-factor obtained 
from the i?-matrix fit of the present data set with literature 
experimental data. The thick continuous line represents the 
extrapolation of the total 5'-factor from the best i?-matrix 
fit, with background poles included and some ANC values 
fixed ("Total i?-matrix"). The dashed-dotted line and the 
thick dashed line are the contributions of the direct capture 
(background poles included) and the high energy resonances 
respectively. All data shown include total uncertainties. 
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TABLE IV. Parameters obtained from the i?-matrix fit 



(MeV) 


(keV) 


p 

i a 

(eV) 


i?^937 


^1121 


^1700 


i ^ \JIiQ V ) 

^2523 ^3791 


^3839 


^4115 


^4652 ^4964 


6.16 (3+) 


14.1 ±0.3 


8± 1 


318 ±4 






47 ± 2 69 ± 2 


140 ±3 


20 ±3 




6.28 (2+) 


11.3 ±0.2 


30 ±4 


930 ± 10 




130 ±15 




175 ±5 






(15) (2-) 


(6 X 10^) 










8.0 X 10^ 








(15) (3-) 


(6 X 10^) 




3.2 X 10^ 












3.4 X 10^ 


(15) (4-) 


(6 X 10^) 




2.6 X 10^ 


2.1 X 10^ 








3.0 X 10^ 





TABLE V. Calculated S-factors 



Ec.m. 


Direct 


Resonant 




Total 




(keV) 


(keV b) 


(keV b) 




;keV b) 




10 


4.92 


0.52 


5.4 ± (th.) 


0.9 ± (exp.) 


0.7 


20 


4.93 


0.54 


5.5 ± (th.) 


0.9 ± (exp.) 


0.7 


40 


4.96 


0.58 


5.5 ± (th.) 


0.9 ± (exp.) 


0.7 


60 


4.99 


0.62 


5.6 ± (th.) 


0.9 ± (exp.) 


0.7 


80 


5.01 


0.68 


5.7 ± (th.) 


1.0 ± (exp.) 


0.7 


100 


5.04 


0.74 


5.8 ± (th.) 


1.0 ± (exp.) 


0.7 


120 


5.07 


0.80 


5.9 ± (th.) 


1.0 ± (exp.) 


0.7 


140 


5.10 


0.88 


6.0 ± (th.) 


1.0 ± (exp.) 


0.7 


160 


5.12 


0.96 


6.1 ± (th.) 


1.0 ± (exp.) 


0.7 


180 


5.15 


1.07 


6.2 ± (th.) 


1.0 ± (exp.) 


0.7 


200 


5.18 


1.18 


6.3 ± (th.) 


1.0 ± (exp.) 


0.8 


220 


5.21 


1.31 


6.5 ± (th.) 


1.0 ± (exp.) 


0.8 


240 


5.23 


1.47 


6.7 ± (th.) 


1.0 ± (exp.) 


0.8 


260 


5.26 


1.66 


6.9 ± (th.) 


1.0 ± (exp.) 


0.8 


280 


5.28 


1.89 


7.1 ± (th.) 


1.0 ± (exp.) 


0.9 


300 


5.31 


2.18 


7.4 ± (th.) 


1.0 ± (exp.) 


0.9 


320 


5.34 


2.53 


7.8 ± (th.) 


1.0 ± (exp.) 


0.9 


340 


5.36 


2.98 


8.3 ± (th.) 


1.0 ± (exp.) 


1.0 


360 


5.39 


3.55 


8.9 ± (th.) 


1.0 ± (exp.) 


1.1 


380 


5.41 


4.32 


9.6 ± (th.) 


1.0 ± (exp.) 


1.2 


400 


5.44 


5.38 


10.7 ± (th. 


) 1.0 ± (exp.) 


1.3 


420 


5.47 


7.00 


12.3 ± (th. 


) 1.0 ± (exp.) 


1.5 


440 


5.49 


9.36 


14.6 ± (th. 


) 1.0 ± (exp.) 1.7 


460 


5.52 


13.24 


18.3 ± (th. 


) 1.0 ± (exp.) 2.2 


480 


5.54 


20.34 


25.2 ± (th. 


) 1.1 ± (exp.) 


3.0 


500 


5.56 


35.65 


40.0 ± (th. 


) 1.1 ± (exp.) 4.8 



In addition, to test the dependence of the extrapolation 
on the adopted fixed ANC values, the latter were varied 
upwards and downwards and the results of the fits were 
compared with each other. For the R/DC^937 transi- 
tion, a 23% variation of the ANC resulted in a change in 
the zero energy extrapolation by 15%. Note that a 23% 
variation of the ANC corresponds to a 50% variation in 
the respective spectroscopic factor. The reported uncer- 
tainties of the spectroscopic factors in [29] are 25%, but 
given the uncertainties involved in converting spectro- 



scopic factors from transfer reactions to ANCs, a conser- 
vative 50% uncertainty was adopted. All the tests men- 
tioned allow for a reasonable upper limit estimation of 
the uncertainty attributed solely to the i?-matrix extrap- 
olation, by adding the variations in quadrature. A 10% 
uncertainty was also added to the extrapolated zero en- 
ergy total S-factor, related to how well the data constrain 
the free parameters. The result for the uncertainty of the 
extrapolation of the total S'-factor to £^ = is ~ 19%, in 
addition to the experimental uncertainty which amounts 
to - 11%. 

In Fig. [lOj the results from this study are compared 
with recent experimental data. The energy dependence 
of the i?-matrix extrapolation agrees very well with the 
low-energy data by Newton et al. [14] (blue inverse trian- 
gles) and Hager et al. [15] (purple stars). In the case of 
the Hager data the agreement with the present extrapola- 
tion is remarkable. The deviation from the Newton data 
is approximately 15% and can be attributed to the uncer- 
tainties of the two experiments (~ 15% for the Newton 
data and 11% for the present data). On the other hand 
the energy dependence of the data from Rolf^ (green 
circles) disagrees with the present extrapolation, given 
that the error bars show the total uncertainty, and an 
overall normalization factor would not change the shape 
of the data. No conclusion can be drawn by the com- 
parison of the present extrapolation and the Chafa et al. 
point (red square) given its large uncertainty. It should 
be emphasized that the literature low-energy data were 
not included in the fit, and thus the extrapolation de- 
pends solely on the present work's data set and on the 
experimental ANCs from the transfer reactions. 

Table [VT] lists the calculated contributions of all the 
measured transitions to the total ^'-factor at zero energy. 
The uncertainties come from the estimated 15% uncer- 
tainty that arises from the choice of the ANC values. The 
experimen tal un certainty is not included. 

In Table VII, the present 5'(0) calculation is compared 
with literature values. The uncertainty of the present 
S'(O) was calculated from the uncertainties of the experi- 



^ The direct capture data by Rolfs are not subject to the normal- 
ization factor that the resonance strengths were corrected by, 
since they were measured relative to -'^^0(p, 7)-'^^F, as discussed 
in section IIIII 
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TABLE VI. S(0) values of ^^0(p,7)^^F for each measured 
7-ray transition 



Transition (keV) 



5^(0) (keV b) 



R/DC^937 

R/DC^1121 

R/DC^1700 

R/DC^2523 

R/DC^3062 

R/DC^3791 

R/DC^3839 

R/DC^4115 

R/DC^4652 

R/DC^4964 



1.7 ±0.3 
0.66 ±0.13 

0.013 ±0.002 
0.17 ±0.03 
0.66 ±0.10 

0.032 ± 0.005 
0.93 ±0.14 
0.55 ±0.08 
0.21 ±0.03 
0.49 ± 0.07 



mental procedure and the uncertainty related to the the- 
oretical fit. The S{0) values for Newton et al. [14 and 
Hager et al. [15 were calculated by fitting the shape of 
the resonant and direct capture contributions as calcu- 
lated from this work to their experiment data, with the 
two parameters being the normalization of each contri- 
bution {Stot = ciiSres-^ci2SDc)- The valuc obtained here 
for the data of Newton et al. (5.1 ±0.9 keV b) is identical 
to the one reported in [14 (5.1 ± 1.1 keV b), but with a 
smaller uncertainty given that part of the previous uncer- 
tainty came from the resonant contribution. No attempt 
was made by Hager et al. [15 to extrapolate the ^-factor 
to ^ = from their data. 

Good agreement is observed with all previous results 
except for the value proposed by Rolfs, which disagrees 
with the present result by a factor of 1.6. The reason 
for this discrepancy is not clear. The evidence given in 
Sec. Ill of this paper for a 40% overall normalization er- 
ror in the Rolfs data (not to be confused with the 0.62 
normalization factor applied to the resonance strengths 
reported by Rolfs [13]) could not solve the issue, as the 
energy dependence of the data is different. For this rea- 
son a definite conclusion cannot be reached. 



TABLE VII. 5'(0) values of ^^0(p,7)^^F from literature 



Source 



S{0) (keV b) 



Rolfs [I3] 
Fox et al. [11] 
Chafa et al. [12] 
Newton et al. [U 
Hager et al. [15] 
Present work 



9.4 ± 1.9 
4.2 ± 1 
6.2 ±3.1 
5.1 ±0.c^ 
5.8±0.c^ 
5.4± (th.)1.0± (exp.)0.6 



^ The listed S(0) value of Fox et al. 11 was calculated by adding 
a resonance contribution of Sres = 0.52 keV b (obtained from 
the present study) to the reported DC value 3.7 keV b. 

^ The Newton and Hager values were estimated by normalizing 
the present resonant and direct capture contributions to the 
respective experimental data. 



The direct capture calculations by Fox et al. [TT] were 
performed with a Wood- Saxon bound state potential and 
hard sphere phase shifts for the scattering wave func- 
tions, and relied on measured spectroscopic factors. The 
result of that calculation was an ^-factor with a slope 
very close to zero. The present extrapolation agrees with 
these results as well, at least for the low-energy region of 
astrophysical interest. The differences can be traced to 
how the calculations were normalized. 

Given the consistency between the experimental data 
by Newton et al.^ Hager et al. [Mj [15], and this work, 
and the relative agreement in the shape of the direct 
capture iS-factor between Fox et al. and this work, it 
is recommended that the direct capture contribution be 
taken as the weighted average of the results from the 
calculated S'dc's from the Newton (4.7 ±0.7 keV b) and 
Hager (5.3 ±0.8 keV b) data and this work (4.9 ±1.1 keV 
hjQ The result of the weighted average is Sdc = 4.8±0.5 
keV b. 



V. REACTION RATES 

The total thermonuclear rate for the ^^0(j9, 7)^^F re- 
action was calculated by direct numerical integration of 
the formula 



Na (crv) = 3.7318 x lO^V 



10,,-l/2rT.-3/2 
^9 



Jo 



(8) 



where the rate is in units of cm^ s~^ mole~^, Tg is the 
stellar temperature in GK, /i is the reduced mass, E is the 
center-of-mass energy in MeV and (j{E) is the reaction 
cross section in barns. However, numerical integration is 
less reliable for narrow resonances where sharp changes of 
the cross section occur. The reaction rate in these cases 
was calculated using the formula for narrow and isolated 
resonances 



Na = 1.540 X 10^\fiTg)-^/\ujj) 

-11.605Er-/T9 



xe 



(9) 



where cu^ and Er are the resonance strength and reso- 
nance energy in the center-of-mass in MeV. The strengths 
for the two low-energy resonances at Er = 65.1 and 183 
keV were taken from Fox et al. [TT] , whereas their energies 
were taken from Chafa et al. [12], following the sugges- 
tion of the reaction rate evaluation by Iliadis et al. [31], 
which is the recommended evaluation by REACLIB [32^ 
for this reaction. The direct capture reaction contribu- 
tion was taken to be the weighted average of the results 
from the analysis of the present experimental data, Hager 
et al.^ and Newton et al.. The direct capture contribution 



Here the theoretical and experimental uncertainties of the direct 
capture part were assumed to add in quadrature. 
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FIG. 11. (Color online) Contributions of different compo- 
nents to the total reaction rate as a function of temperature. 
HE Tails stands for "High Energy Tails" contributions. The 
importance of the direct capture component is evident. 



FIG. 12. Reaction rate ratio of the present reaction rate to 
the REACLIB rate as a function of temperature. The uncer- 
tainty band includes the estimated uncertainties of the direct 
capture and resonant components. 



in the REACLIB evaluation was adopted from the calcu- 
lation by Newton et al. Interference effects between the 
narrow resonances and the direct capture were neglected, 
as the mechanisms mainly proceed through different or- 
bital angular momenta or different 7-ray multipolarities 
than these resonances. 



Figure [IT] shows the percentage contributions of differ- 
ent reaction components to the total reaction rate as a 
function of temperature. "HE Tails" includes the contri- 
butions of the two narrow resonances at Er = 490 and 
531 keV and more importantly the two broad resonances 
at Er = 557 and 677 keV. The importance of the direct 
capture contribution to the total reaction rate is evident. 

The total reaction rates obtained from this work are 



tabulated in Table VIII in the appendix. The ratio of the 
present reaction rate to the REACLIB rate as a function 



of temperature is shown in Fig. 12 The uncertainty band 
includes the uncertainties for the direct capture and the 
high energy resonance tails, as were estimated from this 
work, and the uncertainties of the two narrow resonances 
at Er = 65.1 and 183 keV, as reported in [11 and [12] 
respectively. At the temperature range T = 0.03 — 0.1 
GK relevant for equilibrium hydrogen burning (red gi- 
ants, AGB stars), the reaction rate calculated from this 
work is higher by up to 20%. The reason for this is not 
clear given the fact that the direct capture component 
used by [31] is only by 4% lower than that used for this 
calculation. Nevertheless, there is still reasonable agree- 
ment considering the uncertainties. At the temperatures 
relevant for classical novae (T = 0.1 — 0.4 GK) the present 
reaction rate is also higher than the recommended rate 
from REACLIB by up to 20%. 



VI. SUMMARY AND CONCLUSION 

The excitation functions of the ^^0{p^j)^^F has been 
measured in the energy range of Ep = 365 keV to 1800 



keV at angles ofOiab = 0°, 45°, 90°, and 135° for 10 differ- 
ent primary transitions to states in ^^F. New resonance 
strengths and branching ratios were obtained for the ob- 
served resonances in this energy range. An i?-matrix 
analysis was performed to extrapolate the astrophysical 
S'-factor to lower energies relevant for several astrophys- 
ical scenarios. In this analysis all primary transitions 
observed in this experiment were simultaneously fitted 
together with the ^^0(p, a)^^N data of reference [28] . 
The resulting total S-factor is dominated by the DC con- 
tributions and agrees with recent low-energy measure- 
ments im |T5] within the experimental uncertainties. 

The reaction rate deduced from the i?-matrix analy- 
sis is slightly higher (up to 20%) than the rates in the 
compilations of reference [31 but in agreement within 
the respective uncertainties. The main uncertainty of 
the reaction rate is the direct result of the uncertainty 
of the DC contribution. In the present experimental en- 
ergy range the cross section is dominated by broad and 
strong resonances masking the DC cross sections except 
for the highest energies. The uncertainty could be fur- 
ther reduced by extending the literature data to lower 
energies (< 260 keV). Such measurements are challeng- 
ing and require a low background environment such as 
an underground laboratory. 
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Appendix A: i?-matrix fits performed simultaneously. 



All the remaining i?-matrix fits are shown in Fig. [13 
and 
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Appendix B: Recommended reaction rates 



The reaction rates as a function of temperature calcu- 



Some resonances in energy regions that were 
not covered by the experiment, were also included, when 
they could influence regions where data were taken. All 
fits shown in this appendix and in the main text were shown in Fig. 12 in the main text. 



lated in this work are tabulated in table |VIII| The upper 
and lower rates correspond to the la uncertainty band 
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FIG. 13. (Color online) i?-matrix fits of the differential cross section of ^^0(p,7)^^F. Each plot corresponds to a primary 7-ray 
transition to a different bound state of ^^F. The black squares correspond to 0°, red triangles to 45°, green inverse triangles to 
90°, and blue circles to 135°. The lines are color coded the same as the data points. Only statistical errors are displayed. 
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FIG. 14. (Color online) i?-matrix fit of the ''^0(p, a)^*N Kieser data obtained at 135°. The arrows point to the two 
resonances of interest for the ^^0(p, 7)^*F case. Target effects for this channel were not taken into account. For this reason 
the cross section scale does not apply to the data points near the maxima of the narrow resonances. This simplification has no 
effect on the results for the broad resonances (see text for details). Only statistical errors are displayed. 



TABLE VIII. ^^0(p,7)^*F reaction rates 





Lower Rate 


Reaction Rate 


Upper Rate 






Lower Rate 


Reaction Rate 


Upper Rate 


0.010 


4.14 


X 


10- 


-25 


4.55 


X 


10- 


-25 


4 


97 


X 


10- 


-25 





224 


4.35 


X 


^Q-04 


4.72 


X 


10-04 


5.09 


X 


10-04 


0.011 


6.77 


X 


10- 


-24 


7.45 


X 


10- 


-24 


8 


13 


X 


10- 


-24 





251 


1.13 


X 


^Q-03 


1.22 


X 


10-03 


1.32 


X 


10-03 


0.013 


1.10 


X 


10- 


-22 


1.22 


X 


10- 


-22 


1 


33 


X 


10- 


-22 





282 


2.87 


X 


^Q-03 


3.10 


X 


10-03 


3.33 


X 


10-03 


0.014 


1.43 


X 


10- 


-21 


1.58 


X 


10- 


-21 


1 


72 


X 


10- 


-21 





316 


7.57 


X 


^Q-03 


8.14 


X 


10-03 


8.72 


X 


10-03 


0.016 


1.74 


X 


10- 


-20 


1.91 


X 


10- 


-20 


2 


09 


X 


10- 


-20 





355 


2.22 


X 


^Q-02 


2.40 


X 


10-02 


2.57 


X 


10-02 


0.018 


2.15 


X 


10- 


-19 


2.36 


X 


10- 


-19 


2 


58 


X 


10- 


-19 





398 


7.34 


X 


^Q-02 


8.03 


X 


10-02 


8.71 


X 


10-02 


0.020 


2.30 


X 


10- 


-18 


2.53 


X 


10- 


-18 


2 


75 


X 


10- 


-18 





447 


2.56 


X 


10-°^ 


2.84 


X 


10-°^ 


3.12 


X 


10-°^ 


0.022 


2.23 


X 


10- 


-17 


2.44 


X 


10- 


-17 


2 


65 


X 


10- 


-17 





501 


8.65 


X 


10-°^ 


9.68 


X 


10-°^ 


1.07 


X 


10+00 


0.025 


2.28 


X 


10- 


-16 


2.48 


X 


10- 


-16 


2 


68 


X 


10- 


-16 





562 


2.68 


X 


10+00 


3.01 


X 


10+00 


3.35 


X 


10+00 


0.028 


2.62 


X 


10- 


-15 


2.88 


X 


10- 


-15 


3 


15 


X 


10- 


-15 





631 


7.49 


X 


10+00 


8.44 


X 


10+00 


9.40 


X 


10+00 


0.032 


2.77 


X 


10- 


-14 


3.13 


X 


10- 


-14 


3 


49 


X 


10- 


-14 





708 


1.88 


X 


10+01 


2.12 


X 


10+01 


2.36 


X 


10+01 


0.036 


2.71 


X 


10- 


-13 


3.11 


X 


10- 


-13 


3 


52 


X 


10- 


-13 





794 


4.25 


X 


10+01 


4.80 


X 


10+01 


5.35 


X 


10+01 


0.040 


2.11 


X 


10- 


-12 


2.45 


X 


10- 


-12 


2 


79 


X 


10- 


-12 





891 


8.74 


X 


10+01 


9.88 


X 


10+01 


1.10 


X 


10+02 


0.045 


1.38 


X 


10- 


-11 


1.61 


X 


10- 


-11 


1 


84 


X 


10- 


-11 


1 


000 


1.65 


X 


10+02 


1.86 


X 


10+02 


2.08 


X 


10+02 


0.050 


7.24 


X 


10- 


-11 


8.43 


X 


10- 


-11 


9 


62 


X 


10- 


-11 


1 


122 


2.87 


X 


10+02 


3.24 


X 


10+02 


3.62 


X 


10+02 


0.056 


3.25 


X 


10- 


-10 


3.76 


X 


10- 


-10 


4 


27 


X 


10- 


-10 


1 


259 


4.66 


X 


10+02 


5.26 


X 


10+02 


5.86 


X 


10+02 


0.063 


1.28 


X 


10- 


-09 


1.46 


X 


10- 


-09 


1 


65 


X 


10- 


-09 


1 


413 


7.10 


X 


10+02 


8.01 


X 


10+02 


8.92 


X 


10+02 


0.071 


4.46 


X 


10- 


-09 


5.05 


X 


10- 


-09 


5 


65 


X 


10- 


-09 


1 


585 


1.02 


X 


10+03 


1.15 


X 


10+03 


1.28 


X 


10+03 


0.079 


1.44 


X 


10- 


-08 


1.61 


X 


10- 


-08 


1 


78 


X 


10- 


-08 


1 


778 


1.39 


X 


10+03 


1.57 


X 


10+03 


1.74 


X 


10+03 


0.100 


1.45 


X 


10- 


-07 


1.57 


X 


10- 


-07 


1 


69 


X 


10- 


-07 


1 


995 


1.82 


X 


10+03 


2.05 


X 


10+03 


2.27 


X 


10+03 


0.112 


4.70 


X 


10- 


-07 


5.08 


X 


10- 


-07 


5 


46 


X 


10- 


-07 


2 


239 


2.29 


X 


10+03 


2.56 


X 


10+03 


2.84 


X 


10+03 


0.126 


1.57 


X 


10- 


-06 


1.70 


X 


10- 


-06 


1 


82 


X 


10- 


-06 


2 


512 


2.77 


X 


10+03 


3.10 


X 


10+03 


3.43 


X 


10+03 


0.141 


5.31 


X 


10- 


-06 


5.74 


X 


10- 


-06 


6 


17 


X 


10- 


-06 


2 


818 


3.27 


X 


10+03 


3.65 


X 


10+03 


4.02 


X 


10+03 


0.159 


1.74 


X 


10- 


-05 


1.88 


X 


10- 


-05 


2 


03 


X 


10- 


-05 


3 


162 


3.77 


X 


10+03 


4.19 


X 


10+03 


4.61 


X 


10+03 


0.178 


5.41 


X 


10- 


-05 


5.87 


X 


10- 


-05 


6 


33 


X 


10- 


-05 


3 


548 


4.28 


X 


10+03 


4.74 


X 


10+03 


5.19 


X 


10+03 


0.200 


1.58 


X 


10- 


-04 


1.72 


X 


10- 


-04 


1 


85 


X 


10- 


-04 


3 


981 


4.81 


X 


10+03 


5.30 


X 


10+03 


5.78 


X 


10+03 



